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Abstract: Gravitational tidal forces acting on the virtual e + e~ cloud surrounding 
a photon endow spacetime with a non-trivial refractive index. This has remarkable 
properties unique to gravitational theories including superluminal low-frequency prop- 
agation, in apparent violation of causality, and amplification of the renormalized photon 
field, in apparent violation of unitarity. Using the geometry of null congruences and 
the Penrose limit, we illustrate these phenomena and their resolution by tracing the 
history of a photon as it falls into the near-singularity region of a black hole. 
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1. In general relativity, photons propagate along null geodesies, tracing out the 
causal structure of the background spacetime. Along these null paths, the proper 
time vanishes — in their frame, photons exist only in an instant. They have no history; 
no past, no future, no time in which to evolve or decay. They are caught forever in a 
single instant of time — the unbearable "beingness" of light [1]. 

Quantum field theory, however, profoundly changes this picture. Vacuum polariza- 
tion envelops the photon in a cloud of virtual e + e~ pairs. This dressed photon acquires 
an effective size, given by the Compton wavelength A c = h/mc of the electron, and 
propagates as if it were an extended object, subject to gravitational tidal forces. This 
changes the nature of photon propagation as the phase velocity responds to the gravita- 
tional field and becomes frequency-dependent. In QED, photons see curved spacetime 
as an optical medium with a non-trivial refractive index. 

The propagation of light in curved spacetime therefore exhibits all the usual optical 
phenomena associated with a refractive index n(u), including dispersion and dissipa- 
tion. The size of these effects is 0(a(X c / L) 2 ), where L is a typical curvature scale 
and a is the fine structure constant. However, there are surprises - as first discovered 
by Drummond and Hathrell [2], the low- frequency limit of the phase velocity can be 
superluminal, in apparent violation of causality. Moreover, the imaginary part of the 
refractive index can be negative [3], indicating an amplification rather than dissipation 
of light as it passes through curved spacetime, in apparent violation of unitarity. 

Reconciling these phenomena with the fundamental principles of quantum field 
theory and general relativity has been the focus of our investigations of this subject 
[3-8]. We find that indeed causality and unitarity are respected, but at the expense of 
radical changes to many well-established principles of QFT in flat spacetime, notably 
the optical theorem and the range of theorems associated with analyticity, including 
the vital Kramers-Kronig dispersion relation. In particular, we find that the refractive 
index, now position-dependent, is governed by the whole past trajectory of the photon. 
In QFT, photons are liberated; they have a history, and a rich and fascinating one. 

In this essay, we illustrate all these phenomena by following the history of a photon 
as it falls into a black hole, showing how the superluminal phase velocity is reconciled 
with causality and describing the effects of the gravitational tidal forces on the virtual 
e + e~ cloud. Remarkably, we find that the photon can become undressed — the tidal 
forces squeeze the vacuum polarization cloud, reducing the screening and amplifying 
the renormalized quantum field amplitude as the photon approaches the singularity. 
This, then, is the photon's tale. 
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2. The most illuminating way to describe these effects in QED is in terms of an initial 
value problem for the photon field A^x), which satisfies a wave equation incorporating 
the one- loop vacuum polarization given by the Feynman diagram in Fig. 1. In the 
eikonal approximation, the solution is 



where u is a null coordinate along the classical light ray 7, and V is the corresponding 
null coordinate labelling geodesies in the null congruence around 7. A(uq) is the ampli- 
tude on the initial value surface u = uq, ep}{u) is the polarization vector and riij(u; u) 
is a position, polarization and frequency- dependent refractive index. 



Figure 1. The one-loop Feynman diagram contributing to the vacuum polarization in QED. 

A key observation is that to leading order in (A c / L) 2 , we can replace the background 
spacetime by its Penrose limit [9] associated with the null geodesic 7. We have shown, 
using both worldline and conventional QFT methods [3-6], how the refractive index 
is governed by the geometry of geodesic deviation around 7. But this is precisely the 
property of the spacetime which is encoded in the Penrose limit. This can be seen by 
writing the metric around 7 in null Fermi coordinates and making the Penrose rescaling 
u — y u, v — y X 2 v,x l — > Xx l ; truncating at 0{\ 2 ) then recovers the Brinkmann plane 
wave metric ds 2 = 2dudv — h i j(u)x l x : >du 2 + dx t dx\ with the profile function identified as 
hij = R U iuj [10]. Since the geodesic deviation equation is written in terms of the Jacobi 
fields x % as d 2 x % /du 2 = R l U j U x^ this isolates exactly the required curvature components. 
In QED, the expansion in A is mirrored as an expansion of the Green functions in X c /L, 
which ensures that the gravitational curvature is relatively weak on the quantum scale. 

The refractive index is expressed entirely in terms of geometric quantities related 
to geodesic deviation. Solving the Jacobi equation with "geodesic spray" boundary 
conditions at u = u', as in Fig. 2, determines the Van Vleck-Morette matrix Ajj(w, u') = 
[u' — ^Aj^Ujii'). We have then shown that the refractive index of curved spacetime 
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Figure 2. Null geodesic spray from a point u' on 7. On the left is the equation for geodesic 
deviation of the Jacobi fields x l (u) with a solution written in terms of the 2x2 matrix 
Aij (u, u'). The boundary conditions on the spray are written in the last line. In the illustration 
on the right, geodesies in the spray at u = v! diverge in one direction (blue) and converge in 
the other (red) to cross at a conjugate point. 



in QED with scalar electrons, with initial value surface Uq — > — 00, is 



nij{u;uj) 
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(A similar, more complicated, formula holds for spinor electrons.) Notice the key point 
that riij(u,uj) is given by an integral over the whole of the past worldline u' < u 
of the photon. In curved spacetime, photons have a history, their local propagation 
characteristics being determined by the curvature throughout their path. 

The refractive index has novel analyticity properties, related to the existence of 
conjugate points of the null geodesic congruence. This has far-reaching consequences 
[3] since it shows that amplitudes in curved spacetime have geometric as well as kine- 
matic branch points and cuts, radically modifying many apparently fundamental theo- 
rems in quantum field theory and S-matrix theory. In particular, the Kramers-Kronig 
dispersion relation in curved spacetime becomes simply 



n(u: 00 



, . 1 f°° du . 
n(u; 0) V \ — n[u,cu) . 



In flat spacetime, hermitian analyticity would then imply that the principal value in- 
tegral can be written in terms of Imn(«; u), which would be guaranteed by the optical 
theorem to be positive, forcing n(oo) < n(0). However, causality depends on the wave- 
front velocity (which is identified with the high frequency limit of the phase velocity 
[11]) not exceeding c, that is n(oo) > 1. In flat spacetime, this would clearly be in- 
compatible with a superluminal low- frequency value n(0) < 1. In curved spacetime, 
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however, hermitian analyticity is violated due to the geometric cuts in the complex 
w-plane in n(u; u) and the conventional Kramers- Kronig relation fails. 

Moreover, we also find instances where Imn(«;w) is negative, requiring a reap- 
praisal of the optical theorem itself [7, 8]. The key point is that in curved spacetime, 
the optical theorem, which relates the imaginary part of scattering amplitudes to cross 
sections, only holds globally For example, here we can express the total probability up 
to time u for the decay 7 — >• e + e~ as 

Pit 

P 1 ^ e + e -(u) = Auo I du' lmn(u';u) . 

which is manifestly positive. However, without ^-translation invariance, the local ver- 
sion is no longer necessarily positive and does not represent a true decay rate, removing 
the apparent conflict between lmn(u;u>) < and unitarity. The field may be locally 
amplified, but integrated over its whole past trajectory there must be a net attenuation. 



Imn(u; uj) 

divergence 



\A{u)\ 




u 



divergence 



-> u 



u 



Figure 3. Imn(u;u) and the amplitude |k4(m)| as a function of u in the near-flat spacetime 
region far from the black hole. This shows the characteristic oscillatory transient behaviour as 
the photon becomes dressed with its virtual e + e~ cloud, screening the bare field and reducing 
the amplitude. In QED the initial value of the field is tuned to be divergent to give a finite 
A(u) away from the short-distance region close to uq. 



3. We are now ready to follow the evolution of the photon field as it falls towards the 
singularity of a black hole. Starting from an initial value surface Uq far from the hole 
where spacetime is essentially flat, the imaginary part of the refractive index Imn(n; u>) 
and the field amplitude A(u), for either polarization, have the form shown in Fig. 3. 

Now consider a planar null geodesic in Schwarzschild spacetime with impact pa- 
rameter b. In ref. [6], we developed new techniques for determining the Penrose limit 
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for null geodesies in a general class of black hole spacetimes. We showed that the 
eigenvalues of the plane wave profile function are 

h ± = l -R^k u ±\\^2\^\K S \\ 

where k p is the tangent vector to 7, ^2 = —C py \ p l p m v m p n p = -M/r 3 is the Weyl 
tensor in the Newman- Penrose basis and K s , with \K S \ 2 = 2M _2//3 i) 2 , is the Walker- 
Penrose conserved quantity [12, 13] associated with the null orbits. The Penrose limit 
therefore has hu = —h 22 = 3Mb 2 /r 5 . 

An important simplification arises in the near-singularity limit, with h\\ = — /122 = 
6/25w 2 , independent of M and b. This is a singular homogeneous plane wave. Other 
Penrose limits are readily calculated, e.g. the near-singularity limit of equatorial null 
orbits in the Kerr black hole are also homogeneous plane waves with identical hij, which 
is determined purely by the power-law nature of the singularity [14-16]. 

The refractive index in the near-singularity region now follows from the geodesic 
spray matrix found from the Jacobi equation with this h^, illustrated in Fig. 4. 




Figure 4. The behaviour of the null geodesic spray as the singularity is approached. In one 
plane the geodesies diverge (blue) while in the orthogonal plane geodesies converge with the 
singularity being a conjugate point (red). 

The frequency dependence of n(u; u) for fixed u is shown in Fig. 5. This shows 
explicitly how Ken(u]u) — > 1 as u — > 00 for both polarizations, including the polar- 
ization which is superluminal at low frequency, ensuring causality is maintained. 
The converging polarization also has lmn(u;u) negative, which would violate the 
conventional optical theorem but is consistent with our new curved spacetime formu- 
lation. 
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Figure 5. Ren(u;w) and Imn(ii;w) as a function of frequency ui for fixed position u in the 
near singularity region of a Schwarzschild black hole. The red curves denote the polarization 
perpendicular to the orbital plane, corresponding to a converging null congruence, while the 
blue curves denote the polarization lying in the orbital plane, for which the null congruence 
is diverging. Note that n(it;u>) — > 1 as ui — > oo, as required to maintain causality. 



Figure 6. Imn(u;w) for the two polarizations in QED in a Schwarzschild black hole as u 
approaches the singularity at u = from u < 0. 

We now follow the photon as it approaches the singularity. The amplitude is given 



where we compare with its value at some fixed point U\ in the near-singularity region. 
The results are shown for the two polarizations in Fig. 6. As \u\ approaches zero, we 
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with ci = —0.155 and C2 = 0.105, which implies that the amplitude itself behaves as a 
power law: 

\A {l) (u)\ ~ \u\ ac ^ . 

The fate of the photon is therefore sealed as it falls into the singularity. For 
polarization e^, the curvature pulls virtual e + e~ pairs from the bare field, intensifying 
the virtual cloud and screening the bare field so that its amplitude vanishes. For the 
other polarization, the gravitational tidal forces squeeze the virtual cloud out of 
existence, undressing the photon and revealing the bare field itself. Its history has come 
full circle, from the initial transient dressing and the playing out of renormalization in 
real time as it propagates through curved spacetime, to its final destiny as the vacuum 
polarization cloud once more vanishes as it enters the black hole singularity. 
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